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ABSTRACT 

An  analysis  is  made  of  a  hydrofoil  with  partial  cavitation  in  a 
traveling  gust.  The  problem  is  formulated  in  the  time  domain  with 
pressure  (acceleration  potential)  as  the  fundamental  dependent  variable. 

It  is  shown  that,  in  the  limit,  as  the  aspect  ratio  of  the  foil  becomes 
large,  the  second  time  derivative  of  the  cross-sectional  area  of  the  cavity 
vanishes  regardless  of  the  shape  of  the  gust.  Integra)  equations  are 
presented  through  which  the  length  of  the  cavity  and  the  unsteady  pressure 
distribution  on  the  foil  can  be  determined.  These  equations  are  highly 
nonlinear  and  require  formidable  numerical  procedures  to  sol ve ;  no  attempt 
is  made  to  solve  them  here.  It  is  shown  in  an  appendix  that  the  earliest 
possible  appearance  of  a  nonvanishing  second  time  derivative  of  the 
cross-sectional  area  of  the  cavity  is  to  the  order  of  the  square  of  the 
reciprocal  of  the  aspect  ratio  and  that  to  this  order  the  pressures  at 
distances  large  in  comparison  with  the  chord  of  the  foil  (but  not  the 
span)  grow  logarithmically.  It  may  be  inferred  from  this  result  that  the 
pressures  in  the  far  field  can  be  reduced  by  increasing  the  aspect  ratio 
of  the  foi 1 . 
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NOMENCLATURE 

nondimensional  cross-sectional  area  of  the  cavity 
denotes  wetted  portion  of  foil 
denotes  cavity 

nondimensional  local  chord  of  the  cavity 
c/e 

complex  function  whose  real  part  is  P 

in  the  appendix  H/e  (the  same  as  h  in  the  text) 

£/e 

/T 

unit  vector  in  free  stream  direction 
unit  vector  along  the  chord  of  the  foil 
unit  vector  normal  to  free  stream  direction 
unit  vector  normal  to  _i_ 
a  constant 
a  constant 

nondimensional  local  chord  of  foil 
ordinate  of  cavity  surface 
acceleration  potential 
nondimensional  cavity  pressure 
nondimensional  pressure 
semi -span 

nondimensional  distance  along  span 
time 

nondimensional  time 

perturbation  velocity  due  to  foil 

free  stream  velocity 

nondimensional  perturbation  velocity 

vector  sum  of  free  stream  velocity  and  gust  velocity 

unsteady  velocity  of  fluid 

gust  velocity 

nondimensional  gust  velocity 
denotes  wake 

distance  in  free  stream  direction 
distance  along  chord  of  foil 
nondimensional  distance  along  chord  of  foil 
distance  normal  to  chord  of  foil 


ii 


R-21 18 


I 


» 


\ 


y  nondimensional  distance  normal  to  chord  of  foil 

z  x  +  iy 

a  angle  of  attack 

&  /(A-c)/c 

Y  x/e 

c  small  quantity  inversely  proportional  to  aspect  ratio 

C  mapped  complex  plane  =  £  +  in 

q  imaginary  coordinate  in  mapped  plane 

k  ,<  unknown  functions  of  time 

1  2 

\  dummy  spacelike  variable 

p  nondimensional  difference  between  ambient  pressure  and  cavity  pressure 

v  £/e 

£  nondimensional  dummy  chordwise  coordinate  in  first  section; 

also  in  second  section  real  coordinate  in  mapped  plane 
p  fluid  density 

a  nondimensional  dummy  spanwise  coordinate;  also  in  third  section 

denotes  dummy  time 

r  t/e  ;  also  in  third  section  denotes  t-x 

*  transformed  acceleration  potential 

4>  nondimensional  velocity  potential 

«  y/e 


(  )  except  where  specifically  defined,  denotes  first  order  in  e; 

similarly  for  (  )  ,  (  )3,  (  )4 
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INTRODUCTION 

Ship  propeller  blades  are  known  to  cavitate  intermittently  during 
their  traverse  of  regions  of  large  hull  wake  fractions  or  low  inflow 
velocity.  Such  cavities  are  responsible  for  large  pressure  changes  on 
the  hull  with  concurrent  intense  hull  and  sub-structural  vibration.  It 
is  known  that  the  vibratory  forces  generated  by  cavitating  propellers 
are  much  larger  than  the  forces  generated  by  noncavi tating  propellers. 

The  reason  for  this  is  that  the  cavity  behaves  in  the  far  field  like  a 
source  whose  strength  is  proportional  to  the  second  time  derivative  of 
the  cavity  volume1.  It  is  also  significant  that  cavi ty- i nduced  forces 
arising  at  the  first  few  integer  multiples  of  blade  frequency  are  of  the 
same  order  as  those  at  blade  frequency.  This  is  in  sharp  distinction  to 
the  excitations  induced  by  non-cavi tating  blade  loadings  which  are  in¬ 
variably  small  at  all  higher  multiples  of  blade  frequency,  as  has  been 
pointed  out  by  Breslin1.  Thus,  it  is  clear  that  the  dynamics  of  unsteady 
cavities  on  propeller  blades  must  be  understood  and  reduced  to  a  calcula¬ 
tion  yielding  their  time-varying  geometry  before  the  pressure  field  (which 
induces  the  hull  vibration)  can  be  predicted. 

The  complete  description  of  intermittent  cavitation  on  propeller 
blades  requires  a  three-dimensional  representation  accounting  for  the  mutual 
interaction  of  cavitating  and  non-cavi tating  blade  elements  as  well  as  the 
mutual  interaction  of  the  blades.  This  is  clearly  a  formidable  task.  As 
a  first  step,  we  consider  the  blades  to  be  hydrofoils  acting  independently 
and  we  examine  the  limiting  case  of  large  aspect  ratio.  The  problem  is 
thereby  reduced  to  a  single  two-dimensional  hydrofoil  subjected  to  a 
periodic  gust  pattern  which  is  spatially  stationary  in  hull-fixed  coordi¬ 
nates  and,  therefore,  travels  with  the  negative  speed  of  the  section  in 
foil -fixed  coordinates. 

In  the  past,  there  was,  for  an  extensive  period  of  time,  a  continuing 
controversy  about  the  possibility  of  solving  unsteady  cavity  problems  in 
two  dimensions  because  the  existence  of  time-varying  cavities  of  finite 
length  implies  a  time-varying  source  flow  whose  potential  and  non-convect i ve 
pressure  fields  do  not  vanish  at  large  distance,  but  grow  as  log  r.  Birkhoff2, 
Leehey3  and  others  sought  to  invoke  compressibility  as  a  means  for  resolving 
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this  paradox.  Benjamin4  put  the  question  to  rest  by  pointing  out  that 
such  two-dimensional  problems  are  artifacts  of  attempts  to  reduce  three- 
dimensional  problems  to  manageable  problems  in  two  dimensions  and  that 
they  really  represent  the  inner  behavior  of  three-dimensional  flows  with 
potentials  and  pressures  vanishing  as  a  three-dimensional  source,  i.e., 

/  9  ^  v  “  1/2 

as  (rz  +  s^)  ,  where  s  is  the  spanwise  coordinate.  Thus,  the  logarithmic 

character  of  the  mathematical  representation  of  two-dimensional  time-varying 
cavity  flows  must  be  included  in  order  to  match  the  inner  expansion  of 
the  outer  field.  Benjamin  went  on  to  assert  that  the  rate  of  expansion 
of  the  cavity  section  should  be  regarded  as  a  parameter  of  the  inner  field 
and  could  only  be  determined  by  resorting  to  the  three-dimensional  outer 
field.  With  these  concepts  in  mind,  therefore,  the  problem  will  be  set  up 
as  a  three-dimensional  one  and  the  aspect  ratio  will  then  be  taken  to  be 
large. 

Several  attempts  to  solve  unsteady,  partially  cavitating  hydrofoil 
problems  have  been  made,  but  all  have  fallen  short  of  the  goal  set  here, 
principally  because  they  considered  only  small  disturbances  on  an  estab- 

JU 

lished  steady  cavity."  In  contrast,  we  must  address  flow  patterns  with 
cavities  whose  size  can  grow  and  diminish  over  several  orders  of  magni¬ 
tude  during  one  cycle.  In  one  of  the  earliest  assaults  on  the  non-steadily 
cavitating  hydrofoil,  Steinberg  and  Karp6  studied  the  small  disturbance 
case,  but  they  imposed  a  dubious  condition,  requiring  that  the  velocity  be 
continuous  across  the  foil-cavity  wake  while  ignoring  the  cavity  closure 
condition  altogether.  Wang  and  Wu7  studied  only  small  perturbations  on  an 

q 

established  steady  cavity.  Unruh  and  Bass  presume  that  the  cavity  leading 
edge  can  be  assigned  arbitrarily  downstream  of  the  foil  leading  edge,  whereas, 
many  observations  in  water  tunnels  and  a  vacuum  tank  reveal  that,  on  propel¬ 
ler  blades,  cavitation  almost  invariably  initiates  at  the  leading  edge. 
Furthermore,  they  also  imposed  the,  by  now  familiar,  restriction  to  per¬ 
turbations  of  an  already  existing  steady  cavity. 

The  three-dimensional  problem  posed  by  intermittent  cavitation  on 
ship  propellers  has  been  addressed  by  Johnsson  and  Sdntvedt9,  Noordz i j1 °and , 
more  recently,  by  Kaplan,  Bentson  and  Breslin11.  All  of  these  have  employed 
quasi-steady  strip  theory  utilizing  the  steady,  partially  cavitating,  sec¬ 
tion  theory  of  Guerst1'  with  no  consideration  of  unsteady  flow  mechanics. 

An  exception  to  this  restriction  is  presented  by  Van  Ilouten5  using  a 
numerical  procedure. 
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Despite  this  restriction,  surprisingly  good  correlation  with  measurements 
of  hull  pressures  on  models  and  on  ships  are  reported  in  Reference  11. 

The  problem  is  posed  on  the  pressure  as  the  primary  dependent  vari¬ 
able.  Because  the  length  and  cross-sectional  area  of  the  cavity  (both  un¬ 
known  a  priori)  may  undergo  large  excursions  during  a  cycle,  it  is  not  pos¬ 
sible  to  solve  the  problem  in  the  frequency  domain.  Instead,  it  must  be 
set  up  as  an  initial  value  problem  in  the  time  domain.  As  will  be  seen, 
this  results  in  a  pair  of  nonlinear  integro-d i fferential  equations  (with 
time  as  independent  variable)  which  must  be  solved  numerically.  Inasmuch 
as  the  input  is  periodic,  the  output  must  also  be  periodic  and,  even  though 
an  initial  value  problem  is  posed,  periodicity  will  be  achieved  by  running 
through  several  cycles  of  input  allowing  the  transients  introduced  by  the 
starting  conditions  to  decay.  No  numerical  results  are  presented  here, 
however,  and  the  primary  result  of  the  analysis  is  to  show  that,  in  two 
dimensions,  the  second  time  derivative  of  the  cross-sectional  area  of  the 
cavity  vanishes  identically  regardless  of  the  shape  of  the  gust.  As  a  con¬ 
sequence,  neither  the  potential  nor  the  pressure  grow  logarithmically  at 
large  distance.  In  an  appendix,  it  is  then  shown  that  only  to  the  order 
of  the  square  of  the  reciprocal  of  the  aspect  ratio  is  there  a  possibility 
that  the  second  time  derivative  of  the  cross-sectional  area  does  not  vanish. 
The  solution  at  this  order  allows  for  the  possibility  that  the  three- 
dimensional  flow  behaves  like  a  three-dimensional  source  in  the  far  field. 
The  practical  inference  to  be  drawn  from  this  result  is  that  the  vibratory 
pressures  in  the  far  field  will  be  reduced  as  the  aspect  ratio  is  increased. 
Propeller  blades  of  high  aspect  ratio  are  known  to  be  more  prone  to  cavitate 
than  propeller  blades  of  low  aspect  ratio.  So  a  compromise  must  be  drawn 
between  low  aspect  ratio  blades  which  cavitate  less  and  high  aspect  ratio 
blades  which  induce  lower  vibratory  pressures  in  the  far  field  even  though 
they  are  more  prone  to  cavitate. 

The  authors  are  indebted  to  Program  Director  Mr.  R.  Cooper,  Fluid 
Dynamics  Branch  of  the  Office  of  Naval  Research,  for  support  of  this  work 
under  Contract  N0001 A-77_C-01 22 ,  Task  NR-062-568.  The  analysis  has  been 
conducted  under  Davidson  Laboratory  Project  013. 
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A  PARTIALLY  CAV HATING  FOIL 

Our  present  interest  is  in  the  linear  theory  for  a  partially  cavit- 
ating  planar  foil  which  is  immersed  in  an  inviscid,  incompressible,  fluid 
of  constant  density  p.  We  suppose  that  the  unsteady  velocity  of  the  fluid 
is  given  by 

v  =  v  +  u 

CO 

where  £  denotes  a  perturbation  velocity  while 

v  =  U  i  +  v  (Ut  -  x  )j_ 

— co  —  i  g  i—i 

_i_  =  j_  cosa  +  s  i  not 
1 

j_  =  -j_  s  i  net  +  j_  cosa 
xi  =  x  cosa  +  y  sina 

represents  a  uniform  flow  upon  which  is  superposed  a  traveling  cross  gust 
v  .  We  will  take  v/U  and  the  angle  of  attack  a  to  be  of  the  same  order 
of  small  magnitude. 

Hereafter,  all  terms  will  be  rendered  in  dimensionless  form  by  using 
U;  the  semi-span  S;  and  the  constant  density  p  as  unit  quantities.  The 
linearized  velocity  vector  is  then 

v_  =  _i_  +  [a  +  Vg(t-x)]j+  u_ 

where  j_  and  j_  are  unit  vectors  along  the  x  and  y  axes  respectively.  We 
also  suppose  that  jj  is  source  free  in  the  domain  of  the  fluid  and  that  u_ 
is  irrotational  in  the  fluid  domain  minus  the  wake  so  that  in  the  latter 
domain  there  exists  the  velocity  potential  function  <j>  (x  ,s  ,y ,  t) . 

Figure  1  (see  following  page)  shows  the  area  B  of  the  foil  and  the 
cavity  surface  whose  height  n  is  assumed  to  be  small.  The  trailing  edge 
of  the  cavity  is  assumed  to  lie  always  in  B;  here  x  =  c(s,t)  and  the  foil 
leading  edge  bound  the  area  C  under  the  cavity.  The  quantity  e  which  ap¬ 
pears  in  the  equation  for  the  trailing  edge  of  the  foil  is  supposed  to  be 
small  so  that  the  foil  is  one  of  large  aspect  ratio.  This  foil  with  the 
straight  leading  edge  suffices  for  the  exposition  of  the  basic  ideas. 

Its  wake  W  is  defined  by 
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FIGURE  1.  SKETCH  OF  CAVITATING  FOIL 


W:-  y=0;  l(s)  <  x  f  «;  - 1  <  s  <  1 . 


The  velocity  potent ia 1  <J>  (x , s ,y , t)  must  satisfy 


f>  +  <f>  +  <f>  =0 

xx  Tss  yy 


in  the  domain  of  the  fluid  minus  the  wake  and  be  such  that 


<j>y(x,s,0+,t)  =  -  [a  +  Vg  ( t-x)  ]  c(s,t)<x<H(s);  - 1  <  s  <  1 

<j)^(x,s,0-,t)  =  ”  [a  +  vg  (t-x)  ]  0  <  x  <  SI  (s) ;  -1  <  s  <  1 . 

For  the  same  domain,  a  linearization  of  the  Euler  momentum  equations  pro¬ 
vides  the  first  order  acceleration  potential 

P(x,s,y,t)  *  ♦x+*t--p  +  \*x+*t  +  p\m 

where  p  is  the  pressure.  The  pressure  p[x,s,n(x,s,t) ,t]  at  the  cavity 

wall  is  required  to  be  a  function  of  t  only.  Since  n  is  small,  the  linearized 

version  of  this  condition  is 


p(x,s ,0+, t)  =  pc(t) 
which  leads  to  the  condition 

Px(x,s,0+,t)  +<f>t(x,s,0+,t)  =  -pc ( t)  +  |<{>x  +  4>t  +  pj^  =  \i  (t) 
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for  (x,s)  in  c,  that  is 

0  <  x  <  c(s,t) 

- 1  <  s  <  1 

The  cavity  function  n(x,s,t)  is  determined  by  using  the  linearized  kine¬ 
matic  boundary  condition  for  the  surface 

y  =  n(x,s,t) . 

This  linearized  condition  is 

n  (x,s,t)  +  n  =  a  +  v  (t-x)  +  <J>  (x,s,0+,t),0  <  x  <  c(s,t);  -1  <  s  <  1 
a  t  g  y 

subject  to 

n (0 , s , t)  =0  -1  <  s  <  1 


and  the  closure  condition 

n[c(s,t) ,s,t]  *  0  -1  <  s  <  1  . 

Let  us  begin  the  determination  of  <(>(x,s,y,t)  with  the  representation 


<j>(x,s,y,t)  =  t—  //  - - - 777  *  [<j>]d£da 

^  B+W  [  (5-x)2  +  (ct-s)z  +  y2]3/2 


-7 ~  //  - - - TTV  ’  [<f>„]dCda 

B+W  [  (£-x)2  +  (0-s)2  +  y2]  y 


t«J>]  =  <i>(£,a,0+,t)  "  <(>  (C  »o  »0-  ,  t) 

[4>y]  =  4>y  (C  ,cr,°+,  t;)  -  $  (c,o,0-,t) 

This  implies  a  representation  for  the  acceleration  potential  P(x,s,y,t). 
We  have 


(—  +  -~)  | 

' 3x  3t;  1 


"T - //  - * - 5/7*  t<^>] dCda 

^  B+W  H  [(C-x)2  +  (a-s)2  +  y2] 3/2 


IT  //  - - - 077  •  t4> d^da 

^  B+W  [  (C-x) 2  +  (a-s)2  +y2]3/2  1 
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3/2 

At  the  leading  edge  [<J>]  =  0;  therefore,  since  y/[  ]  -*■  0  as  £  -*■  «, 

an  integration  by  parts  yields 


_  1 
1  Air 


// - * - 

B+W  [ (C-x) 2  +  (a-s)2  +  y2] 


372 


(  E^]  +  E<t>t]  )d£da  . 


However,  we  require  the  acceleration  potential 
P(x,s,y,t)  =  <(>x  +  4>  t 
to  be  continuous  in  the  wake  and  so 


(JL  +_L)  |  =  _Lrr 

'3x  3t''l  knJL 


E(C-x)2  +  (o-s)2  +  y2] 


[P]dCdo  . 


Since  we  also  require  the  normal  velocity  $  to  be  continuous  in 
the  wake,  we  have  [<}>^]  =  0  there,  from  which 


'2  =  -iW/ 


1 


^  B  [ U-x)2  +  (a-s)2  +  y2l  ^2  y 


[<j>  ]d£da 


1  «.(a) 

if  I 


.  I  (o) 

TT  /  4  ] dx  d£da 


- 1  o 


f  (C~x) 2  +  (a-s)2  +  y2] ' /2  H  t-  y 


=  i L 

l»ir  * 


£(a) 


- 1  [x2  +  (a-s) 2  +  y2]  o 


172  /  ^ Jd5  da 


,  1  £(a) 

if  I 


(g-x) 


-1  o  [(c-x)2+(a-s)2  +y2] 


*-(a) 

-  3/2  *  /  (Od*  d5d<* 

+  v213/Z  £  y 


Then 


£(a) 


-I  [x2  +  (a-s)2  +  y2] 


772  *  /  t*y]dC  d° 


i  *  *(<*)  a  r 

1  /  / 


’"-l  o  85\[(E-x)! 

1  (a) 


(g-x)  ] 

+  (a-s) 2  +  y2] 3/2J 


+  i/  / 


(g-x) 


ih) 

I  [<f>  IdA  d£da 

g  Y 

H(a) 

■£r  f  [*  ]dX  dCda 


-1  o  [(£-x)2  +  (a-s)2+y2]3/2  8t  5  y 
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and  an  integration  by  parts  of  the  second  term  on  the  right-hand  side  of 
the  last  equation  shows 


(-1.  +  JL) , 

' 9x  9  t  1  2 


-  J_  JL  [ 

4tt  9t 


1 


_ 1 _ 

[x2  +  (a-s)2+y2]1/2 


i(a) 

/  Cd>  I dC  do 

A  ' 


1  l(a) 

/  / 


-1  o 


_ (g~x) _ 

[(g-x)2  +  (a-s)2  +  y2] 3/2 


[<{.y]dx- t<J>y] 


ji-dCda 


Downstream  of  the  partial  cavity  and  on  B,  we  have  [4>y]  =  0  by  virtue  of 
the  boundary  conditions.  Therefore 


<£  +  £>'2 


J_  _9_  }  _ 1 _ 

3t  -1  [x2  +  (o-s)2+y2],/2 


c(cr.t) 

J  [41  ]d£-do 

A  ' 


1 


I 

o  [(C-x)2  +  (a-s)2+y2]3/2l3t 


(g-x) 


f  9 

U-  f  u 


]dX-[4>  ] 

y  y 


-} 


d£da 


This  can  be  expressed  in  terms  of  n(x,s,t)  if  we  use  the  condition 
nx(x,s,t)+nt  =  <j>y(x,s,0+,t)  +  cx  +  v^(t-x)  *  [4>y]  which  holds  on  C;  and  the 
closure  condition 


n[c(s,t) ,s,t]  =  0 

wi  th 

n  (0  ,s ,  t)  =  0 
We  find 

c(o.t)  c  .  c 

/  U  ]dg  =  /  ntd?  =  —  /  n(g,o,t)dC  *At(a,t) 
o  y  o  o 


where  A  is  the  cross-sectional  area  of  the  cavity.  Also 


a  c(o,t) 

M(g,c,t)  H  —  /  [4>y]dX  -  [4>y] 


gY  /  {n^ (X ,0 , t)  +  n{}dX  -  n^(c,o,t)  -  nt 


L 
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as 


3  c(°*t) 

M(f,,o,t)  *—  /  nt(X,o,t)dX  -  2nt(fa,o,t)  -  n^(C,o,t)  . 

Our  development  shows  that  the  acceleration  potential  can  be  expressed 


P(x,s,y,t)  =  ♦  (x.s.y.t)  +  * t (x , s  ,y , t) 


l*ir  ■' 


-1  [x2  +  (o-s)2  +  y2] 


172  *  Att 


A  (o,t)do 


+  c(V}  (e-x)  M(s.o,t)dedg 

^-1  o  [  (c-x)2  +  (o-s)2  +  y2]  *  2 


1  l(a) 

il  I- 


y [P]d^da 


1  o  [  (£-x)2  +  (o-s)2  +  y2] 


3/2 


which  is  the  same  as 


P(x,s,y,t)  «  -57  Att^S,t)  *n(*2+y2) 


+  7—  •—  j  S-L  g-n [  |o-s)+  /(0-5)2  +  x2  +  y2]A  (o ,  t)do 
Air  3  s  a  -  s  1  tt 

+  _LJ_  /  C(°/t)  (C-x)  (s-o)  M(C.o.t)dCdo 

4,r  3s  -1  o  (£-x)2  +  y2  [(s-o)2  +  (C-x)2  +  y2] 1/2 


Ltt  A*.  1  j 


(s-o)  [P]dF,do 


4ir  35  -I  o  U-x)2+y2  [  (s-a)2  +  (^-x)2  +  y2]  1/2 

We  introduce  the  quantities: 

£  =  ev  ;  c(o,t)  *  ec(o,t) 


x  =  ey 

t  *  ET 


;  A  -  e2  A 


y  =  cm  a,  (0)  =  eh(o) 

I 

which  implies  a  high  aspect  ratio  wing  and,  because  of  the  scaling  of  the 
time,  it  also  implies  oscillations  whose  wave  length  is  of  the  order  of 
the  chord. 

We  now  seek  an  asymptotic  expansion  for  4>(x,s,y,t)  such  that 
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I 


5 


<fr(x,s,y,t)  -  $ (cy ,s ,eu),eT) 

where 

,  IhillfL  o 

Le-*o  v^lc) 


I  e\(cHk(  Y.s.w.t) 
k=  1 


Each  is  to  be  determined  by  equating  terms  of  like  order  after  the  ex¬ 
pansion  has  been  substituted  for  $  in  the  representation  for  P.  However, 
for  the  time  being,  let  us  forego  a  detailed  development  of  the  expansion 
and  be  concerned  in  this  paper  with  only  the  lowest  order  term  4>, (x  ,s  ,y ,  t) 
In  the  appendix,  however,  the  higher  order  terms  are  examined. 

Apart  from  the  term 

T^Att(s,t)  in  (x2  +  y2 )  ,  (  =  £  in(x2+y2)), 

the  lowest  order  approximation  to  P  can  be  obtained  by  neglecting  the  com¬ 
binations  x2  +  y2  and  (£-x)2+y2  which  appear  in  the  radicals  in  the  last 
expression  for  P.  With  respect  to  the  non- integrated  logarithmic  term, 
we  need  to  consider 


"5tT  Att  tn  (y2+w2) 

and  hence 

7=  Altt(*.t)  “/  [♦ e  . 

o 

Without  loss  of  generality,  we  can  take  vj(^)  *  1  so  that  the  order  of 

27  Aitt(s’t}  in  e 

is  given  by  in  e.  Now  this  order  cannot  be  matched  with  that  of  any  other 
term  which  appears  in  the  expansion  of  P.  We,  therefore,  conclude  that  we 
are  forced  to  set 

Altt(s,t)  "  °  -I  <  s  <  1  . 

With  this,  it  follows  that  the  first  order  approximation  to  P  is  given  by 
Pj  (x,s,y,t)  -  $u(x,s,y,t)  +  $lt(x,s,y,t) 
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,  3  1  c(o^t)  s_o  (C-xjMjtS.o.tldCdo 

Uv  3 S  J  ■»  I  <?-rrl  "  TT  - 


-1  O 


(C-x)2+y2 


i  ^  1  (cr) 

ii  f  s-g 

Air  3s  ,  •*  1  s-g 

- 1  o  * 


yfPjJdedo 
(C-x)2  +  y2 


,  c(s.t)  (C-x)M  (C,s,t)d€ 
P  (x,s,y,t)  =  —  / - 1 - 

o  (C-x)2+y2 


,  Us)  y  [PjJdC 

2it  o  (C-x)2  +  y2 


We  see  from  the  last  result  that,  with  s  as  a  parameter,  Pj(x,s,y,t)  is 
a  two-dimensional  harmonic  function  in  the  exterior  of  the  cut  along  the 
x-axis  from  0  to  £(s);  and  we  also  see  that  Pj  vanishes  as  /xz  +  y2  -+■  «>. 

Analytic  Determination  of  ^ 

Let  us  turn  to  the  determination  of  Pj  =(f>lx(x,s,y,t)  +  $u  and  ^(x.s.y.t). 
For  typographical  convenience,  we  will  drop  the  subscript  and  often  suppress 
the  parameter  s.  With  this  understanding,  the  first  problem  we  have  to 
solve  is  this:  Find  P(x,y,t)  such  that 

P(x,y,t)  =  <f>x(x,y,t)  +  <f>t  (x  ,y ,  t) 

is  harmonic  in  the  exterior  of  the  slit  along  the  x-axis  from  0  to  fi. ( s)  , 
while 

P(x,0+,t)  =  4>x(x,0+,t)  +  <j>t(x,0+„t)  =  p(t)  0  <  x  <  c(s,t); 

Py(x,0+,t)  =  4>yx(x,0+,t) +(jiyt(x,0+,t)  =  0  c (s , t)  <  X  <  i(s) ; 

Py(x,0-,t)  -  4yx(x,0-,t)  +<j>yt(x,0-,t)  =  0  0  <  x  <  l(s)  . 

We  assume  that,  if  P(x,y,t)  possesses  any  singularities,  then  they 
exist  at  only  one  or  both  of  the  points 

(0,0)  ;  (c,0+)  ; 

where  (0,0)denotes  the  leading  edge  of  the  foil  and  (c,0+)  denotes  the 
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moving  trailing  edge  of  the  cavity.  Let  us  suppose  that  the  physically 
acceptable  solution  for  P{x,y,t)  is  the  one  with  the  mildest  singularities 
at  (0,0)  and  (c,0+).  With  respect  to  the  behavior  in  the  neighborhood 
of  infinity,  we  impose  the  condition  dictated  by  the  three-dimensional 
analysis,  namely 

P  -*■  0 

as  /xz  +  y^  00 

If  z  =  x  +  iy,  the  function  defined  by 


maps  the  exterior  of  the  slit  along  the  x-axis  from  0  to  1  into  the  quarter 
plane  Im  5  >  0,  Re  £  >  0  where  Im  denotes  the  imaginary  part  and  Re  denotes 
the  real  part.  For  this  mapping,  the  branch  /(z-z) /z  is  defined  as  shown: 


/{Z~z)/z  =  |  /(Z~z)/zj  exp  i(01-0)/2 


y 


FIGURE  2.  PHYSICAL  PLANE 
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(The  quantity  6  =  /(£-c)/c) ,  for  example,  is  negative.)  Then  Lj  is 


mapped 


M 


2 


M 


3 

O- 


FIGURE  3-  MAPPED  PLANE 

into  Mx  defined  by 

£  *  0  ;  «.  >  n  >  o  ; 

Lj  is  mapped  into  M2  defined  by 
n  *  0  ;  0  <  £  <  /^5  ; 

and  L  is  mapped  into  M  defined  by 
n  -  0  ;  /^6  <  E  <  oo 

The  point  at  infinity  in  the  z-plane  is  mapped  into 


J l  — 


-0  -  /P7 


Under  this  mapping  P(x,y,t)  is  transformed  into  a  function  4>(C,n)  which 
is  harmonic  in  £  >  0;  n  >  0  with  singularities  at  C  -  0  and  c  =  ®.  The 
function 

4>  -  p(t) 

vanishes  along  Mj  while 


3n 


-  u(t) } 


vanishes  along  M2  and  M3  .  By  using  the  theory  of  analytic  continuation, 


13 


R-21 18 


it  follows  that  $  is  the  real  part  of  a  function  of  5  which,  except  for 
poles,  is  analytic  in  the  S-plane.  (The  continuation  shows  that  we  cannot 
admit  branch  point  singularities.)  If  we  admit  only  the  mildest  singular¬ 
ities  at  Z  =  0  and  C  =  “,  we  have 


4>  =  u { t)  +  ^(ORe  ~  +  *2(t)Re  Z 
or  in  terms  of  z 


P(x,y,t)  =  y  (t)  +Kj(t)Re 


1 


l-z 


+  <2(t)R eyy-^2--  3 


-  6 


In  order  to  avoid  excessive  writing  in  the  sequel,  let  us  express  the  last 
result  as 


P(x,y,t)  =  Re  F (z , t) 

We  now  have  the  problem  of  finding  the  velocity  potential  <j>(x,y,t) 
which  is  harmonic  in  the  xy-plane  minus  the  positive  x-axis;  exhibits  a 
wake  downstream  along  the  x-axis  from  x  =  i  to  x  =  and  which  satisfies 

<f>x  +  t  =  Re  F(z*0 

The  general  solution  of  this  equation  can  be  expressed  as 


o 

<t>(x, y,t)  =  /  Re  F(o+z,  o+t)do  +  4>(x-t+a,y,a) 

a-t 

where  a  is  an  arbitrary  constant.  We  may  suppose  that  <j>(x,y,0)  =  0  for 
all  (x,y).  Then,  if  we  take  a  =  0,  we  have 

o 

4>(x,y,t)  =  /  Re  F(a+z,a+t)d0 
-t 

This  function  is  harmonic  for  0  <  ang  z  <  2ir;  and,  if  we  write  it  in  the 
form 

x 

<t>(x»y,t)  *  /  Re  F (X+iy,X+t-x)dX  , 

x-t 

it  is  not  difficult  to  verify  that  it  shows  a  downstream  wake  in  the  sense 
that  <(>  and  ^suffer  a  jump,  provided  t  is  sufficiently  large,  as  z  crosses 
the  downstream  x-axis. 

lit 
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Determinat ion  of  the  Coefficients 

The  next  step  is  to  enforce  the  two  boundary  conditions  which  4> 
must  satisfy:  the  closure  condition,  and  the  condition  on  the  regular 
part  of  P  at  infinity. 

The  component  $  is 


g 

f>  (x.y.t)  =  /  r-  Re  F(cr+z,c+t)do 

y  _t 


=—  /  Re  i  F(o+z,a+t)do 
aX  -t 

a  1 

=  — /  Re  i  F (a+z-t ,o)da 
o 


For  0<x-0i<i.,we  must  have 


[a+v  ( t-x) ]  «  <j>  (x,0-,t)  =  —  /  Re  i  F(a-t+x-0i  ,o)do 

g  y  dx  Q 


t-x  t 


=  —  {  /  +  /  Re  i  F (a-t+x-Oi ,c)do} 

dX 


o  t-x 


However,  Re  i  F (a-t+x-Oi ,a)=  0  for  0  <  a  -  t+x  -  Oi  <  l  or  t-x  <  a  <  t  +  £-x 
end,  consequently, we  must  have 


3  t_x 

-[o+v  (t-x)]  =  g —  /  Re  i  F(a-t+x-0i ,a)dcr 
9  o 

Then,  if  we  set  t-x  =  x ,  we  see  that  for  t  >  0 


a  +  vg(x)  =  - —  {  Re  i  F(cr-x,a)da 


or 


T  T 

/  [a  +  v  (a)]da  *  /  Re  i  F(a-x,a)da 
o  9  o 


(1) 
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This  is  the  first  equation  which  relates  the  unknowns  k  ,  k  and  8. 

1  2 

For  c(t)  <  x  +  Oi  <  we  must  have 


~[a+v  (t-x)]  =  /  Re  i  F (o-t+x+Oi  ,o)da 

o 

a  t-x  t 

=  -^  {  /  +/  Re  i  F (o-t+x+Oi ,o)do} 

X  o  t-x 

We  already  know  that 

3  t-x 

-[a+v  (t-x)]  =  t—  /  Re  i  F(a-t+x-0i  ,a)dcr 
X  o 

3  t_x 

=  —  /  Re  i  F(o-t+x+0i ,a)da 
o 

for  t  -  x  >  0  because  the  space  variable  o-t+x,  in  the  range 
- (t-x)  <  a  -  t+x<0 

required  in  the  last  integrals  lies  ahead  of  the  foil  and  F  is  continuous 
there.  Therefore 

a  1 

0  =  —  /  Re  i  F(a-t+x+0i ,o)do  (2) 

t-x 

c  ( t)  <  X  <  l 

is  the  second  equation  which  relates  the  unknowns. 

Equation  (2)  can  be  simplified.  A  change  in  the  dummy  variable  of 
integration,  a-t+x  =  X,  gives 


Re  i  F (X+Oi ,A+t-x)d^ 


c  <  x  <  l 


(2a) 


Note  that 

Re  i  F(X+0i ,X+t-x)  =  0 

when 


c(X+t-x)  <  X  <  l 
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This  inequality  defines  an  upper  limit  for  the  integral  in  (2a).  For 
simplicity,  suppose  that 

A  =  c(A+t-x) 

has  but  one  root  A^  for  each  t-x.  Then 
A  =  c(X1+t-x) 

and  this  root  is  a  function  of  t-x  so  that  we  can  write 
Xj  =  A^t-x) 

(The  root  Ai  is  given  by  the  ordinate  of  the  point  of  intersection  of 

X  =  c(a) 
and  the  1 i ne 

X  =  a  -  t  +x 

which  passes  through  the  point  (t,x)  with  slope  equal  to  I.)  Consequently, 
equation  (2a)  becomes 

Xx (t-x) 

■—  f  Re  i  F(X+0i ,X+t-x)dX  =  0 
X  o 

or 

c (Xj+t-x) 

■—  /  Re  i  F (X+Oi , A+t-x)dA  =  0  ; 


and,  if  we  introduce  t-x  =  T,  we  have 
c(Xj+t) 

4-  !  Re  i  F (X+0 i ,X+T)dX  =  0 

dx 

o 

Integration  gives 
c(Xj+t) 

/  Re  i  F (X+0 i , X+x)dX 
Xj(0) 

=  /  Re  i  F(x+0i,X)dA  =  cons't.  =  K 

o 
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This  can  be  written  as 
c(Xj+t) 

/  Re  i  F[A+01  ,A+t+A  -c(A.+-r)]dA  =  K 
o  1  1  ° 

and  if  we  set 

t"  =  Aj  +  t 

we  have 

c(t“)  * 

/  Re  i  F[A+0i,A+t  -c(t")]dA  =  K  (2b) 

o 

A 

where  t  is  arbitrary. 

We  now  turn  our  attent  ion  to  the  cav  i  ty  surface ,  the  equation  for  which  is 
y  =  n(x,s,t)  0  f  x  ^  c(s,t) 

or,  if  we  suppress  s, 

y  =  n(x,t)  0  <  x  f  c(t) 

and,  as  we  have  seen  in  the  three-dimensional  analysis,  n(x,t)  can  be  deter¬ 
mined  by  using  the  linearized  kinematic  condition 

nx (x '  t)  +  nt(x,t)  =  4>y  (x,0+,t)  +  [ci+vg  (t-x)]  0  <  x  <  c(t) 

In  terms  of  F  this  is 

3  t 

nx(x,t)  +  n^ (x , t)  /  Re  i  F(a-t+x+0i ,o)da  +  [a+vg(t-x)] 

3  t-x  t  0 

=  —  /  +  /  Re  i  F (a-t+x+Oi ,a)do}  +  [a+v  (t-x)] 

X  o  t-x  9 

which,  by  virtue  of  the  first  condition  (1),  reduces  to 

3  t 

n  (x,t)  +n.  (x,t)  =—  /  Re  i  F (a-t+x+Oi  ,o)da  0  <  x  <  c(t) 

X  t  oX  .  “ 

t-x 

for  t-x>0.  Now  n(x,t)  must  vanish  at  the  leading  edge,  i.e., 
n  (0, t)  =  0 

and  the  only  solution  of  the  differential  equation  for  n(x,t)  which  satisfies 
this  condition  is 
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n(x,t)  =  /(x- A) Re  i  F (A+0 i , A+t-x)dX 

o 

The  formula  for  the  cross-sectional  area  A(s,t)  determined  by  the 
cavity  surface  and  the  foil  has  already  been  used  in  the  three-dimensional 
analysis.  Here  it  is: 

c(t) 

A(t)  =  /  n(x,t)dx 

o 

c(t) 

=  /  [c(t)  -  A]Re  i  F[X+Oi ,A+t-c(t)]dA 

o 

The  closure  condition  to  be  imposed  is 
n  [c  ( t)  ,  t]  =  0 

This  can  be  enforced  by  integrating 


nx (x , t)  +nt(x,t)  =—■  /  Re  i  F(a-t+x+Oi  ,a)do 


t-x 


=  —  /  Re  i  F (X+Oi , X+t-x)dX 
3x  o 


0  <  x  <  c ( t) 


Since  n(0,t)  =  0,  we  have 


n(x,t)  +  /  nt(c,t)d£;  =  /  Re  i  F(X+0i  ,X+t-x)dX 


or,  using  n[c(t),t]  =  0, 


c(t) 


c  ( t ) 


c(t) 


/  ~  n(C,t)dC  =  /  n(C,t)dC=  '/  Re  i  F (X+Oi  ,X+t-c)dX 


3 1  dt 

o  o 


which  is  the  same  as 


c(t) 


A ^ ( t )  =  /  Re  i  F[X+0i, X+t-c(t)]dX 


(3) 
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This  equation  relates  the  rate  of  change  of  the  cross-sectional  area  of 
the  cavity  to  the  unknown  parameters  k  ,  and  8.  However,  Eqs.  (2b) 
and  (3)  imply  that 


from  which 


(as  we  already  know  is  required)  and 

A  =  K  t  +  K 
o 

If  a  monotonic  growth  of  the  cavity  cross-sectional  area  is  precluded, 
it  follows  that  we  must  take 

K  =  0 
o 

It  appears  then  that  the  first  order  approximation  produces  a  constant 
value  for  the  cross-sectional  area  of  the  partial  cavity.  Equation  (3) 
then  becomes 


c(t) 

/  Re  i  FU+Oi  ,X+t-c(t)]dX 


o 


0 


(3a) 


and  we  now  have  equations  (1)  and  (3a)  for  the  determination  of  ,  <2 , 
and  6. 

The  last  equation  for  the  determination  of  the  unknowns  is 


L  [  zl-x-  ReF  (2’t}  =  0 


or 


Kl(t)  _ 

Re{p(t)  +  +  k  (t) /i  -  8)  =  0  . 

/t^b  2 


(4) 


The  equations  (l),  (3a)  and  (4)  form  a  set  for  the  determination  of 
kc i ,  <2  and  8.  The  equations  are  linear  in  Kj  and  <2,  but  disconcertingly 
nonlinear  in  8-  Two  of  them  are  highly  nonlinear  integral  equations  in  t 
and  the  parameter  s.  Their  solution  for  given  a,Vg(t)  and  y(t)  can  only 
be  accomplished  by  using  formidable  numerical  procedures. 
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The  result  4>x  +  <f>t  =  Re  F(z,t)  presented  on  page  l^*,  where  F(z,t)  is  f 

also  given  there,  can  be  shown  to  subsume  some  classical  results.  For 
example,  if  the  dependence  on  time  is  eliminated,  and  Kj  and  <2  are  eval¬ 
uated  from  the  kinematic  condition  that  the  vertical  velocity  is  simply 
a,  then,  together  with  the  (steady)  closure  condition,  the  result  can  be  § 

shown  to  reduce  to  that  of  Guerst12. 

As  another  example,  take  Kj  =  0  and  write 

1/2 

*x+*t-v(t>  +  v^^Re  (1  -  —  /z)  f 

As  £3-*-°°,  take  p  -  /-§"  k2  =  0  and  we  have 

i 

<2 

4>x  +  4>t  =  —  Re/U-z)  /z 

which  is  the  case  of  the  unsteady,  completely  wetted,  foil. 

The  analysis  presented  thus  far  is  in  a  somewhat  unsatisfactory 
state  inasmuch  as  it  has  been  shown  that  A^t  =  0  so  that  no  two-dimensional 
source-like  flow  at  distances  large  in  comparison  with  the  chord  of  the 
foil  is  possible.  It  is  clear  that  this  result  comes  about  because  the 
aspect  ratio  has  been  taken  to  be  infinite.  It  is  desirable  to  determine 
to  what  order  in  aspect  ratio  a  non-zero  value  of  A  will  occur.  To  this 
end,  we  are  devising  a  procedure  for  the  methodical  determination  of 
higher  order  asymptotic  approx imat ions .  If  e  is  the  reciprocal  of  the 
aspect  ratio,  and,  if  the  orders  of  the  asymptotic  approximations  are  * 

characterized  by  1,  e,  t2Z ne,  e2....,  we  find  that  there  are  strong 
indications  that  A  remains  equal  to  zero  at  least  until  k  corresponds 
to  c2.  Thus,  it  may  be  expected  that,  as  the  aspect  ratio  is  increased, 
the  time  variations  of  the  cavity  cross-sectional  area  and,  hence,  the  • 

variation  of  the  total  volume  of  the  cavity  on  the  three-dimensional 
foi 1 )  will  diminish. 

Finally,  comparison  of  the  formalism  of  this  report  with  the  two- 
dimensional  analysis  of  Van  Houten  shows  that,  whereas  we  have  demonstrated 
the  necessity  for  taking  A  =  0,  thereby  eliminating  the  source-like 
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behavior  of  the  pressure  at  large  distances.  Van  Houten  has  retained 
this  term  in  his  work.  We  believe  this  to  be  improper,  especially  in 
view  of  the  requirement  that  the  boundary  condition  at  the  cavity  can 
only  be  a  statement  that  the  difference  between  the  cavity  pressure  and 
the  ambient  pressure  is  given,  whereas,  if  the  pressure  grows  logarith¬ 
mically  at  large  distances,  no  ambient  pressure  can  be  defined. 

The  solution  of  the  integral  equations  which  Van  Houten  analyzes 
numerically  is  equivalent  to  finding  the  acceleration  potential 

P"(x,s,y,t)  =  +  <t>'t 


=  17  Att(s,t)  + 


*n2|a_s!  \t(s’t)da 

.  ,  c(y}  u-x)M*dc 

+  27  J  Tc:xTr+y'7 


where 


,  I  £(rs)  y[P"]d€ 

+  I7  /  Tf-xFTy 

A  c(s.t)  * 

Aj.  (s , t)  =  /  [<f>  jd£ 


M" 


3t 


c(s^t) 

5 


JU 

[<jT]dX  -  [4>y] 


and  P  is  subject  to  the  cavity  condition 

P  (x,s,0+,t)  =  u(t)  0  <  x  <  c(s,t);  -1  <s<l; 
and  the  condition 

P"(x,s,0±,t)  =  0 

along  the  wetted  surface  if  we  take 
4>y(x,s  ,0^,t)  =  -a  -  vg{t-x) 
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there.  The  representation  for  P  shows  that  it  may  possess  a  logarithmic 
singularity  at  infinity.  However,  if  we  admit  such  a  singularity  when 

JU 

P  is  determined  by  the  function  theory  method  used  above,  it  turns  out 

JU 

again  that  Att(s,t)  must  be  zero  in  order  to  satisfy  the  boundary  condition 
on  the  upper  side  of  the  foil  downstream  of  the  cavity  and  the  closure 

JU 

condition.  Thus,  it  appears  that  the  first  order  approximation  A  =  0 
can  be  regarded  as  a  consequence  of  the  imposed  boundary  conditions  for 
the  two-dimensional  approximating  problem  rather  than  a  consequence 
of  matching  in  an  asymptotic  expansion. 
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We  proceed  to  give  a  method  for  the  theoretical  determination  of 
higher  order  approximations  for  the  cavitating  foil.  It  is  beyond  the 
scope  of  the  present  report  to  supply  all  of  the  details  of  the  develop¬ 
ment.  It  is  expected  that  a  fuller  analysis  will  be  produced  later.  Here, 
our  object  is  to  find  the  order  of  approximation  which  gives  the  first 
indication  that  Att  nowhere  A  is  the  cross-sectional  area  of  the  cavity. 

We  have  found  that  the  acceleration  potential  is 


P(x,s,y,t)  =  <(ix(x,s,y,t)  +  4>^(x,s,y,t) 


=  _  J_  (±  +  ±)  } 


[<{>^]d£da 


kir  '*x  '  3t'  -1  o  [ (£-x)2  +  (a-s)2  +  y2] 1/2 

+  J_  }  £(ro)  _ y[Pld£do _ 

^  1  o  [  (S~x)2  +  (a-s)2  + y2] 


where 


[<t>]  =  [♦(C.o.t)]  =  ♦(?,a,0+,t)  -  <J>(C,a,0-,t) 


In  order  to  satisfy  the  boundary  condition  for  the  cavity,  we  must 
satisfy 

1  1  9  9  l  c(V)  [<Ud£da 

p(x,s,o+,t)"Y  [P(x,s,t)]  =  -  j—  (gr+Tr)  /  f  - - - “7/2 

2  Ait  3x  3t  o  [  (c-x)2  +  (o-s)2]  1/2 


for  0  <  x  <  c(s,t)  ;  -1  <  s  <  1 

In  order  to  satisy  the  boundary  condition  for  the  wetted  part  of 
the  foil,  we  must  satisfy 


=  +  I  (— +  — ) 

~  2  v3x  3t' 


[<j>y(x,s,t)] 


1  32  }  l(?> 

^4TT  3x3S  , 


/(x-g)2  +  (s-a)2 [P]d£da 
(x-C) (s-o) 


An  integration  gives 


±  j  $v(£,s,0±,t)dC  +  ~  /  4»v(C,s,0±,t)dC 


3t 
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7  {^T  +  7$  I  I*/*-*, t)]dc 

X  1 


+  _LJ_  f  Mf0)  {SU<E  *  (s-g)2 

4  dS  ’  -  (x-O(s-o) 


- 1  o 


- }  [P]d£do 

s-o 


and  if  we  solve  this  for 


/  <j>  U,s,0±,t)dE  we  find 
x  y 


/  *  (5,s,0±,t)dS-±  4- J  t<f>  (C,s,t)]d£ 

v  '  w  ' 


+  i  I  h  /  /  ^^(x^Tls-of^  '  fp(C,a,X+t-x)]dCdcrdA 


4tt  j  3  s  . 
x-t  -1  o 


In  dimensional  terms,  the  normal  velocity  along  the  wetted  part  of 
the  foil  is  required  to  be 

-  [lla  +  vg  (Ut-x)  ] 

In  dimensionless  terms,  this  velocity  is 


-[a  +  -f  (St-Sx)l 


where  S  is  the  semi-span.  For  a  foil  of  large  aspect  ratio,  we  can  take 


S 

e 


where  e  is  small  and  k  is  the  chord  length,  so  that  the  normal  velocity  can 
be  expressed  as 


or 


-(«  <t-x>]> 


-[« *  ;,<*£)] 


In  dimensional  terms,  the  acceleration  potential  at  the  cavity  wall  is 
P(t)  . 

A-2 
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In  dimensionless  terms,  it  is 


><t>  ,  p<§> .  - , 

PUZ  PUz  P'eJ 


In  dimensional  terms,  the  trailing  edge  of  the  foil  is 
x  =  Sit(s)  =  Sef(s) 
and  in  dimensionless  terms 
x  =  ef (s) 

In  dimensional  terms,  the  trailing  edge  of  the  cavity  is  given  by 
x  =  c ( s , t )  =  Sec (s  ,  -Sj~)  =  Sc(s,t) 
and  in  dimensionless  terms 

x  «  c  c(s,  jjy)  =  ec(s,~) 

Since  the  unknowns  depend  on  the  parameter,  e  ,  let  us  introduce 
x  *  ty  ;  K  =  ew  ;  y  =  ey  ;  t  =  et 
and  assume  the  inner  asymptotic  expansions 

oo 

4>(x,s,y,t)  =  0  (ey ,s ,ey  ,ct)  =  e$(v,s,y,t;e)  =  [  cv^e)  1*>IC (v .s  ,u  ,t) 

K=1 


c(s,~)  *  l  vk(c)c<(s,t) 

e  k*1 


where 


v  +  1 

L  — — rx-\-  m  0 

e-*-o  v  (e) 


In  terms  of  the  new  variables,  the  boundary  condition  for  the  cavity 
£(v,s,0+,t ;e)  -  y  fP(y,s,T;e)] 

,  .  ,  1  c(o,t)  (w,a,i;e)]dci)do 

*  -U—  f  f  _ H - 

^  ^  3T  -1  o  [  (s-a)2  +  r2] 1/2 

A-3 


where  r2  =  g2(w-y)2 


and  0  <  y  <  c 

This  can  be  expanded  into 


P(Y,s,0+,T;e)  -  j  [P(y>s ,t ;e) ] 


=  (-L  +  J-) 


c(s.t) 


—  /  Anc  j  aj~Y  |  ‘ 


,  „  1  c(o,x)  I  I 

7 k&l  I 


-1 


s|  •  [$  ]duda 

y 


(A“ 


<\j 

c 


“  8tT  Ts*  /  e2(w"Y)2^n  e|u-y|  ‘  UJdu. 

o 


J_  33  i  CKt) 

Utt  3 sJ  A  ■*  / 

-1  o  < 


a-s 


a-s 


-(u-y)2 


-  '  -Tp(w-Y)2^n2|o-s|  j.  |f ^Jduda 

+  o(e2;a,s,w,Y) 


In  terms  of  the  new  variables-  the  boundary  condition  for  the  wetted 
surface  is 


'Xj 

.  c 

1  t  r'V 


/  <f>  (w,s,0±,T;e)du  =  ±  j  f  [<p  (w,s  ,t; c) ]du 

Y 


,  1  }  —  f  f  ^  e/(s-a)zVq2  £  r r>  / 

L-rr  '  3<;  v  '  {  q  (s-o)  "*  s-n  J  *  * 


4tt  j  3  S  . 

Y"T  -1  O 


s-a 


,X+t-y;c) ]dwdadX 


where 


q  =  c  (X-u) . 

An  expansion  of  this  yields 


°°  'v  ]  C  -v 

/  <(>  (aj,s,0+,x;e)dw  +  -r  /  [<J>  (to,s  ,t  ;e)  ]dn 

Y  P  Y  u 


c  }  »  }  f(rg)  [i*(u,.  o ,  X  +t  ~y  '>  c )  ]  dwdod  X 


R-2118 


_  _j_  y  [P  (<ij  ,s  ,X+t~y;£)  ] dcodX 

?TT  1  >  X-W 


y"T  O 


-  3  c2 

4  s -a 


.  Y  »3  1  f(tf) 

/  /  -f 

Y-T  -1  o 


r  s-a 


s-a 


£n2  ]  s-a  | 


■e2  h  (a-s)  2-ne  |  X-d 


y  (X-w) [P (w , s ,X+t-y; e) ] dwdadX 

(A— 2) 


1  yt  33  -  f(?} 


+  TJf  /  /  /  o(e2  ;a ,  s ,X ,w)  •  [?  (w.a.X+x-y;  e)  ]dwdadX 

Y-t  -1  o 


where 


h(a-s)  = 


a  >  s 


a  <  s 


The  expression  (A-2)  shows  that  the  inner  asymptotic  expansion  of 
4>(x,s,y,t)  should  have  the  form  <j>(x,s,y,t)  =  <j>  (ey  ,s  ,ey  ,ex)  =  e  $  (y,s  ,p  ,t  ;e) 
where 


£(y,s,u,x,e)  =  <1>1(y,s,vj,t)  +  e  <P2  +  e2ine 

Then  for  P(x,s,y,t),  we  have 


>3  +  g24>4  + 


P(x,s,y,t)  =  4>x(x,s,y,t)  +  <t>  t  (x  ,s  ,y ,  t)  or 

S'  *\j  .  -x,  .  . 

P  =  4> yC  Y» S  ,V  ,t;e)  +  <f>T  ( Y,s ,p ,t  ;e) 


=  P  (y,s,p,t)  +eP2 


+  e22.ne  •  P  +  e2P,  +  . . . 

3  *♦ 

The  function  P(x,s,y,t)  is  a  three-dimension  potential  function  in 

^2  ^2  ^2 

the  exterior  of  the  foil.  If  we  apply  the  operator  (-r— +  — 7-  +  — 7-)  to  the 

ox  o  o  y 

assumed  asymptotic  expansion,  we  find 


A-5 
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l  v  (e)  [P  +  P  +  c2P  ]  =  0 
<  <yy  i<yy  <ss 

and,  hence,  each  P^(Y,s,y,x)  must  satisfy 

P  (v,s,u,t)  +  P  =  H  (y,s,p,t) 
kYY  xyp  k 

in  the  exterior  of  the  slit  defined  by 

0  <  y  <  f(s)  ;  v  =  0 
For  example,  we  have 

P  +P  =0  ic  =  1,2,3 

kyy  Kyy 

while 


P,  +  P,  =  -P, 

4YY  ‘♦UP  lss 

The  equations  for  the  4>  's  can  now  be  found  by  equating  terms  of 
like  order  in  (A-l)  and  (A-2) . 

Consider  first  <f>  ^ :  In  (A-l)  there  is  no  match  for  the  term 
Ene  ,  3  3  Cl(s/x)  . 

~  (37 +  j  [ny]du  ; 

o 

and,  therefore,  we  must  take 
c. 


+  =  &  +  T&  Aix<s»T>  =  Awt(s»t)  =  0  • 

o 


•l 


3y  3x  {  v3y  '  3x'  "Itt’ 


Here  A  is  the  first  order  approximation  to  the  cavity  area. 

1 

For  the  cavity,  4> ^  must  then  satisfy  the  integral  equation 

Cj (s,t) 

p(t)  -  J  [P^Y.s.t)]  =  27^  +  ^  /  An | cu— Y |  !♦,  ]«!«  i 

o  u 

and,  for  the  wetted  surface,  it  must  satisfy 

•  1  C)  ,  x  f(s)  [P  ((o,s,X)]du)dX 

/  4>  (w,s  ,0±  ,  x)  doi  =±  —  J  [ 45  (u,s,x)]dw  +  ■=—  /  /  - 

Y  IP  2  Y  lw  2  o  o  X +  Y-  T  -  u) 

from  which 
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.  i  a  T  f(s)  [P,  (w;s,X)]dwdX 

(y.s.o+.t)  -  *  }[*1„(v.s.,)!  i  ! 

o  o 


^lUY^’5’01’7^  +  <t>1yT^Y,S,0±,'t^ 

1  s  ~  f(s)  [P, (w,S,T)]dw 

=  ±2(w+37)Uiy(Y’s>T)3 '^rw  l  ~~r~^ — 

where,  on  the  wetted  surface, 

<t>  1^  =  '[a  +  v  (t-y)] 

As  we  have  seen,  these  integral  equations  can  be  solved  by  first 
finding  the  analytic  function  F(z,t)  such  that 

Re  F  (z ,  t)  =  <(>  lx  (x ,  s  ,y ,  t)  +  4>1(;(x,s,y,t) 

,  a  a  C1  , _  _  ,  f(s)  wtPjdu) 

=  27(w+~}  /"n/(“  y)  +u"  [\]du  +  27  /  T — ^“TI 

o  ^  o  (w- y)  +  P 


l  a  a  c(s,t)  _____  ,  *(* 

r>  /  h.jdt+i  / 


,  Us)  y[Pj]d? 
ly^d^+2ir  I"  (C_x)z  +  y^ 


1  c  (s  ,  t )  (c-x)  Mi(c,s,t)dC  1  Us)  y[Pj]dC 
2tt  ^  (C-xU  +  yz  2tt  ^  (%-x)2  +  yz 


The  above  representations  and  the  integral  equations  imply  that  we  must 


Re  F  (x+0 i , t)  =  p ( t) 


for  the  cav i ty  and 


=  0 

y=0± 

<t>ly(x,s,0±,t)  =  -  [a+Vg  ( t_x)  ] 

along  the  wetted  surface  while  Re  F  (z,t)  vanishes  at  infinity. 

It  has  been  shown  that  the  first  approximation  to  the  cavity  surface 
y  =  n  (x,s,t)  can  be  found  from 
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nlx(x,s,t)  +  nlt(x,s,t)  =  [4>  (x ,  s  ,0,  t)  ] 


=  /  Re  i  F(X+0i,  X+t-x)dX 


=  Re  i  F (x+0 i ,t)  +  /  ~  Re  i  F(X+0i ,X+t-x)dX 
o 


for  0  <  x  <  c.  From  this 


c-e 


[<t>  (c-e,s,0,t)]  =  Re  i  F(c-e+0i , t)  -  / 

o 


~  Re  i  F(X+0i , X+t-x) 


x=c-e 


Re  i  F(c-e+0i  ,t)  1  9  C  rE  „  :  _  _  .  , 

-  -  “  - t—  /  Rei  F(X+0i  ,X+t-c+e)dX 


.  •  ,  .  9t 

1  -  c  1-c  o 


if  e  is  a  sufficiently  small  positive  quantity. 
The  closure  condition 

n (c , s  ,  t)  =  0 


g  i  ves 


/  n  ^  t  ( £,  s ,  t)d£  =  /  [4>ly(C,s,0,t)]d^ 


iy 


or 


A  =  / Rei  F(X+0i  ,  X+t-c)dX 


or 

3  C 

A  itt  =  /Rei  F (  X+Oi  ,  X+t~c)d  X  ; 
o 

but  we  have  found  that  A  must  vanish.  Therefore 


‘W')>ly(c'£’S’0’t)]  =  Lc^o 


Rei  F  (c-e+0  i  ,  t) 
1  -  c 


which  is  infinite  as  can  be  seen  from  the  explicit  expression  for  F.  On 
the  other  hand 


Lr^o  Uj  (c+e,s,0,t]  =  0 
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Consider  second  <ji  :  For  the  determination  of  <J>2 ,  we  need  to 
consider  the  term 


2tt 


£ 


c  /  [.}>  ]dw 


We  leave  c  unexpanded  and  define  it  to  be  equal  to  c^.  Indeed, 
[$  ]  vanishes  beyond  so  that  such  terms  can  make  no  contribution 

to  <J>2 •  A  similar  argument  will  hold  for  the  higher  orders.  Since  no 
other  term  of  order  e  fi-n  e  appears  in  (A— 1 ) ,  we  must  take 


'■'l 

A„  =  -r-  /  U,  ]  do  =  0 
2tt  3t  j  2p 
o 


We  see  now  that,  even  for  the  second  order  approximation,  the  second 
derivative  of  the  cross-sectional  cavity  area  is  zero. 


With  this,  <f>2  must  satisfy 


P,(y,s,0+,t)  -  y  [p2(y,s,t)J  =  -57  (^7  +  -^r)  /  ftn|u>-y|  [<)>,„] du> 

o 


c^s,  ) 


2'” . .  2  1 2v,’J'l,J  “  2tt  v3y  ’  St7  i  "  LT2pJ 


for  the  cavity. 

For  the  wetted  surface,  <t>2  must  satisfy 
c, 

oo  1 

/  <j)  (u,s,0±,x)du  +-  /  [<(|  (u,s,T)]du> 

Y  M  y 


x  1  f(o)  [P.(w,a,X)]do)dadX  t  f(s)  [P, (w,s,X) )dmdX 
I  f  o  f  f  1 


k  TT  '  dS  .  J 

o  - 1  o 


s-a 


+  J_  /  / 

2*  i  o  *+Y-t-u) 


from  wh i ch 


_  1 


(>2(j(y»s,o±,t)  ;  j  [^2p (y»s*t)  1  = 


,  a  t  f(s)  [P  (co,s,X)]dwdX 

-  J-  -If  [  — i - 

2tt  3y  J  J  X+y-T-w 
o  o 


and  hence 


t>2yY  (y ,  s  ,0+  ,t)  +  <J>2ijt(y,s,0±,t) 


,  a  a  la  f(,s)  [p,  (u,s,T)]dw 

o 
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These  integral  equations  can  be  solved  by  first  finding  the  potential 


funct ion 


Re  F2(z,t)  =  <(>2x(x,s,y,t)  +<J>2 1  (x ,  s  ,y ,  t) 


=  J_  (—+—) 

2w  v3y  3r; 


,  f(s)  y[P2]dw 

^2y]du  +  27  J  T 0,-VF  +  7 

o 


±  (-L+-L) 

2 it  '3x  3t; 


_ _ _ _  ,  y [P2]dC 

-x)  +  y  [4>2y]^  +  27  [  Tr-xT7"+"7 


c  (C-X)M,  d?  ,  l(s)  y[pjd? 


„  l  [  ■  z  -  _i_  ■ 

2tr  ■*  (?-x)z  +  yz  2 it  * 

o  o 


-x)z  +yz 


The  function  Re  F2(z,t)  must  vanish  at  infinity  and  the  integral 
equations  show  that  it  must  be  such  that 

Re  F2  (x+0 i , t)  =  P2(x,s,0+,t) 

for  0  <  x  <  c,  while  along  the  wetted  surface 

~  Re  F2(z,t)  [  =  <t>2xy (x ,s ,0, t)  +  ♦  (x.s.O.t) 


where  (x ,s  ,0 ,  t)  is  an  induced  velocity  on  the  wetted  surface  and 
P2(x,s,0+,t)  is  an  induced  acceleration  potential  at  the  cavity.  The 
latter  quantities  are  to  be  determined  by  matching  inner  and  outer  expan- 


Consider  next  4>  3 :  The  component  3  is  to  be  determ  i  ned  by  equat  ing  terms 
of  order  e2  2,n  e  in  (A- 1 )  and  (A  -  2 ) ;  but  first  we  need  to  consider  the  terms 

C1 

■s—  (~  +  -r~)  I  [4>  ]dw 

2  n  3y  3t  ^ 

Because  the  term  of  order  e2(ilne)2  cannot  be  matched  in  (A-l),  it 
follows  that  we  must  have 


a  -..a-/, 


S,,  "3,  ■ 0  • 
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We  can  now  see  that  $  must  satisfy 
P 3 ( Y , s , 0+ , t )  -  y  [P 3(y»s, t) ] 

r 


,  1 

—  /  *n|o>-Y| 

O 

i  Cl 

\  +  2?  [  [4>4y]d“ 


U) 


1  Jt 

157  3s': 


(to— y ) 2  [4>ly]dw 


for  the  cavity,  i.e.,  0  <  y  <  c , ;  y  =  0+. 

For  the  remainder  of  the  foil  <j>3  must  satisfy 


/  <i>  (a),s,0±,T)dw 

Y  M 

,  C1  ,  T  f(s) 

=  ±  7  /  [4> 3  (u,s,T)]do)  +Yjj-  /  /  [P3(w,s 

Y  o  o 


.  T  „3  1  f(o) 

'  tst  /  ir?  /  /  h(o_s)  *  tpi(u 

o  - 1  o 

from  which 

,  1  a  T  f(s)  tp3(w 

(Y ’ s »0± *t)  -  ±  2  [♦3w(y.s,t)1  -  27  37  /  /  “177 


t  3  1  f(a) 

+  77  ^  Ts7  /  ^  h(°-sHP1(^.~;?)]rf-dadA 

o  - 1  o 

'*>3^Y,S,0±,T^  +<^3pT^Y»s»0±,x) 

-!T<W*W>  C*3u(y,s.t)] 


A)  ]daidA 


a, A) ]dudadX 


s ,A) ]dudA 
x-u) 
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J_  _3_ 

2n  3y 


Us) 

/ 


[Pg  (<1>,S  ,A)]dw 
Y"o) 


3  1  f (?) 

+  4*  TT*  f  f  h(o-s)  [P1  (o),a,X)]dudo 
-Jo 


These  integral  equations,  which  express  the  boundary  conditions, 
can  be  solved  by  first  finding  the  potential  function 

P3(y,s,p,t)  =  Re  F3(z,t) 

C1 

=  2*  (w  +  f  £nA“-v Y4+u*  [4>3p]d“ 

O 

]  f(s)  p[P3]dco  ,  g 

+  2i r  £  (w-y)^  +  +  2ir  3t  &  ^4p]du) 

-sr  <w  +  £>i&  /  «“-*>2+“2}  [♦,„]«» 

o 

1  32  f(,s) 

”47  3^  /  wfPjIdiu 

o 

where  z  =  Y  +  ip.  From  the  integral  equations,  we  find  that  Re  F3(z,t) 
must  satisfy 

Re  F3(y+0i,T)  =  P3(y,s,0+,t)  for  the  cavity  and 


37  Re  F3{z't) 


p=0± 


=  <f>3lJY(Y,s,0,T)  +  <^>3ijt(y  ,s  ,0,t) 


for  the  remainder  of  the  foil;  while  F,(z,t)  must  have  a  pole  at  infinity. 


If 


J_  _3_ 

2ir  3x 


/  [(tl4y]du 


is  not  zero,  <J>4  remains  to  be  determined  from  the  equations  for  order  e2 
and  there  is  a  coupling  between  the  approximations  of  third  and  fourth 
order.  In  the  solution  for  ReFg,  the  induced  quantities  ^^(YiS.O.x) 
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and  P3(y,s,0+,t)  are  at  first  unknown.  Ultimately  they  are  to  be  found 
by  matching  inner  and  outer  expansions. 


Finally,  let  us  consider  4>u :  Along  the  cavity  interval  <j>  must 


satisfy 


P4(y,s,0+,t)  -  j  [P4(y,s,t)3 


r  ci 


27  /  an | ai“Y |  t^4  (u, s ,t) ] da) 


■  11  f  I  o-s 

/ITT  V  J 


=  (_L  +  _L)  / 

v3y  3t;  ' 


1 

/ 

-1 

c. 


an  2  j  o-s  |  /  [<t>Ull]  doidCT 

o 


Air  3s  o-s  i  lT4y 


"87^7  /  ^"Y)2  an ] o)“Y |  Uj^du) 


(w-y)20  +  2  |  o-s  I  } 

- 1  o 


[<(>1^]duda 


Along  the  wetted  surface,  <t>4  must  satisfy 


(y>s  ,0±,t) 


.  l  a  T  f(s)  [Pu(u,s,X)]da)dX 

-  t  4  K„hr.»,T)l  -££/  / 

o  o 


U7J  3s  . 
o  - 1  o 


from  which 


7  s-a  2  s-a 

+  h(o-s)  { £n |  X+y-t-o)|  +1 } 


\Pj  (w,o,X)]da)dodX 


^i...v^Y*s,0±,t)  +  iJ'4ijt{y.s,0±,t) 


4py 


T  1+.., Cy.s.t)]  I 


1  3 


f(s)  [Pu(io,s,T)]du 


2  v3y  3t  iy4P 


27r  3y 


Y~a> 


+  4-  A-  )  f(/a)  f  3lsi2l+  >  hzsl 

kv  3sT  jj  o  (  *  s-a  ^  s-o 


an  2 | s-a ] 
+h(o-s)  {  fcnjywl  +  1} 
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> [Pj (w,a,T)]dwda 


In  accordance  with  these  integral  equations,  the  function 


Re  FJz.t)  =  (^7+  -^)  /  In/fo-Y)*  +  U2  [♦4Jd“  +  97 

o 


'l 


I 


4  '*■ ' 1 '  2 IT  '3y  '  3t '  i  ■-  ir4 yJ“”  2 IT 


'1 


2ir 


_ _  a  ,  f(s)  y [P.ldto 

tn/Yz  V  yr  /  [*»„ld“  +  2V  / 


dw 


must- satisfy 

ReF4(Y+0i,T)  = 

P  (y.s.O+.x) 
4 


/  ^■q~ s-l-  £n  2 1 a-s |  •  /  [<f>  ]dwda 
-1  a-s  o  “u 


_1 _ 3_  r  ia-s 

Air  3s 


+  {W+^T)V^^/  (w-y)2  Hn | ui  y |  [♦jjdu, 


1  C1 

+  17-t4  /  (“-Y)2  (1+tn  2 1 a-s |  }  [* 

,  loir  9sJ  \  '  o-s 

V,  "•  0 


for  the  cavity,  and 


3lTRe  F4(Z'T) 


y  =0± 


=  *4PY(y’S,0±’t)  +  <J>4yT(Y,S,0±,T) 

-X&  }  '*/*(! iaz-L *  i  iazaL ta  2|« 

An  3s3  1  4  s-a  2  s-a 

o 

|+  h(a-s)  {£n  |  y-4u  |  •**  1) 
for  the  wetted  part  of  the  foil. 


[F 


f(s)  ytPiJdw 
£  (u>-y)z  +  U2 


,  Jdwdcr 
ly 


^dtoda 


As  we  have  seen,  the  function  P4 


6,  +  must  satisfy  the  Poisson 

ay  4t 


equat i on 


P 

4yy 


+  p. 


4UP 


P 


lss 


where 


P^Y.s.u.t) 


-L  (  J-t±) 

2tt  V  3y  3t;  J 

o 


£n /Ty-to) 2  +  y2  t^^ldu) 


,  f(s)  y  f  Px  3  dto 
2tt  q  (y-u) 2  +  y2 


It  can  be  shown  that  if  is  an  appropriate  particular  solution  of  the 
Poisson  equation,  then 

P4  =  S4  +  Re  F4(z«t) 

The  potential  function  Re  possesses  a  logarithmic  singularity  at 
infinity,  namely, 


271)7  /  [<\y]du)  *  Zn/^  +  y2 

o 


'I 


provided 


a  1 

A  =  —  /  [<{>  ]du 
4tt  3t  j  4y 


does  not  vanish.  So  far  there  is  no  apparent  reason  for  supposing  that 

A,  =  0. 

4TT 

We  conclude  from  our  analysis  that,  if  A(s,t)  is  the  cross-sectional 
cavity  area  and  if  A  does  not  vanish,  then  it  must  be  of  order  equal  to 
or  higher  than  e2. 
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